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Abstract 

The contribution of interference effects to the radiation of a charged particle moving in a medium of 
randomly spaced plates is considered. In the angular dependent radiation intensity a peak appears at 
angles 9 ~ ir — 7 _1 , where 7 is the Lorentz factor of the charged particle. 



I. INTRODUCTION 

It is well known that a charged particle passing through a stack of randomly spaced plates radiates 
electromagnetic waves (see, for example, [1]). The radiation is caused by the scattering of the electromag- 
netic field (pseudo-photon) of the charged particle from the inhomogeneities in the dielectric constant. 
In an earlier study one of us has shown [2] that, in analogy with three-dimensional random media [3], 
the multiple scattering of the electromagnetic field plays an important role. However, in the multiple 
scattering approach only the diffusion contribution was taken into account. At this level the approach is 
equivalent to the radiative transfer theory for light transport in e.g. slab geometries, see [4] for a recent 
review. 

On the other hand, interference effects are important when waves propagate in random inhomoge- 
neous media. Anderson localization [5] and the enhanced backscattering peak [6] are manifestations of 
these effects. Other interference effects show up in correlations and higher moments of the transmitted 
intensity. They were also reviewed in [4]. 

In the present paper we want to investigate interference effects for radiation of a charged particle 
moving in a system with one-dimensional randomness. 

II. FORMULATION OF THE PROBLEM 



The system which we want to study consists of a stack of plates randomly spaced in a homogeneous 
medium. The dielectric constant of the system can be represented in the form 

e(z, to) = e (u) + J2 I 6 H - e oMl Il e (^ a / 2 ) - ®( z ~ z i + «/ 2 )l] > W 

i 

where Zi are the random coordinates of the plates, a is their thickness, and e (co) and b(ui) are dielectric 
permeabilities of the homogeneous medium and the plate, respectively. It is convenient to represent the 
dielectric permeability as the sum of an average and a fluctuating part 

e(z, lj) = e + e r (z, w), < e r (z, w) >= 0, (2) 

where e —< s(z,oj) > and averaging over the random coordinates of plates is determined as follows 

<f(z,uj)>=[l[pf(z,z i ,w), (3) 

i z 

where L z is the system size in the z-direction. The vector potential of the electromagnetic field created 
by a moving charged particle satisfies the equations. 
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divA e(r, u)ip{r, uj) = 0, 

c 

uj 2 

W 2 A + —e(f,uj)A(f,uj) = j(r,uj) (4) 
cr 

where j(f, uj) is the current associated with the moving charged particle 

j(r,uj) = -— -S(x)S(y)e^/\v\\z (5) 
c v 

The symmetry of the problem allows us to choose the vector potential along the z-axis, Ai = SziA. The 
electric field is expressed through the potentials in the following way 

E(r,uj) = — A(r : uj) — grad<^(r, uj) (6) 

As usually, we decompose electric field into two parts: E = Eq + E r . Here Eq is the field originated 
by the charged particle moving in the homogeneous medium with dielectric constant e and E r is the 
radiation field associated with the fluctuations of dielectric constant. The radiation tensor is determined 
as follows 

hj(R) — E ri (R)E* 3 (R) (7) 

where R is the radius- vector of the observation point which is far away from the system, i? S> L. For 
expressing the radiation intensity through the radiation potential A r , we decompose the vector potential 
analogous to the decomposition of electric field A = Ao + A r . The fields Aq and A r satisfy the equations 

c 2 ' 

2 2 9 
Vi, + -oSA r + -^ErAr = je r A (8) 



V 2 A + ^jeA =j(r,oj) 



Using (Q) and (|6|) one can express the radiation tensor in terms of radiation potential 

< I,j(R) >= < A r (R,u,)Al(R,u,) > + 5 -S- < A r (R.cj)-g^L-A;(R,u) > 

For obtaining (Q) we supposed that e r <C e. In order to carry out the averaging over the random 
coordinates of the plates it is convenient to express the radiation potential in terms of the Green's 
function of equation (||) 

J 2 

r-2 



A r (R) = ~ I e r (r)A (r)G(R\r)dr 



c 
uj 2 

V 2 + k 2 + -irEriz) 



G{r,r') = 5{r-f l ) (10) 



where k = ujy/e/c. 
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III. GREEN'S FUNCTION 



In this chapter we follow the approach of our previous paper [2]. For the bare Green's function in 
momentum representation one has from ( |To| ) 

In the coordinate representation one has 

^^-^-/^ (12) 
The average Green's function in the independent scatterer approximation has the form 

° {g) = fc2-y» + »ImS(g) ' (13) 
where the imaginary part of the self-energy ImE(g) is determined self-consistently by the Ward identity 

dp 1 



( 27r ) 3 " u/fc 2 - ql 



B(i?.-V* a -9pi)+ fl (i?*+v* 3- ^ |) j' Wpl<k (14) 

Here B(p) = (2ir) 2 5(p p )B(\p z \), where p" p is the transverse component of the p and B{\p z \) is the 
correlation function of the one-dimensional random field 

B{\z-z'\) = ^ <e r {z)e r {z')> (15) 

Using explicit form (]]]) of the random field and carrying out the averaging with the help of (g), we obtain 

4(b-e) 2 nsm 2 q z a/2 ^ 
= ~2 74 ( 16 ) 

where n is the concentration of plates in the system. The photon mean free path in the z direction is 
determined by 

'» - W" (17) 

As one could expect, the photon mean free path depends on the momentum direction. This is principally 
different from the isotropic case. In the case when the momentum is directed along z, using ( jl3 ) and 
(|l7|), one obtains 

<< 9 = °> = swfW*) (18) 



3 



From this point we shall call this quantity the pseudo-photon mean free path. It follows from ( fuf ) 
that B(Q) = (b — e) 2 Lo A na 2 /c 4 and B(2k)/B(0) ~ sin 2 ak/a 2 k 2 . Therefore we have for the mean free 
path 

f 4k 2 /B(0), fca> 1 
l = l(9 = 0)& { ' W (19) 
[ 2k 2 /B(Q), fca< 1 

Note that the expressions (|l9|) are obtained in the Born approximation which holds when | ^Jbje— 1 1 ka <C 
1. Notice also the restriction imposed on the angles. As was mentioned in [2] our consideration is correct 
up to angles 9 « vr/2 - 8 (where \8\ > (l/fc0 1/3 )- 



IV. THE RADIATION INTENSITY 

Substituting ( |l0| ) into (|^) we obtain following expression for the averaged radiation tensor 

< Iij(R) >= J drdf'A (r)A*(r') < e r (r)e r (r>) 

> (20) 



8 M 8^G(R,r)G*(r',R) + " " " " + ^^G*--^- + 

c b crs^ oR l oz oRjdz e c 4 oR l oz e c 4 oRjdz 

As mentioned above, the observation point is far away from the radiating system. For this reason, using 
([l2]), one can obtain following useful relations 

Gn(R r) m L e <fc(«-«r) d 2 G {R,r) k 2 mn z ik(R _ m 

G °^' rj ^i? 6 ' dR t dz 4nR 6 '^ >>r W 

Here n is the unit vector in the direction of the observation point R. The radiation tensor (|2(]) consists of 
three contributions. Single scattering and diffusion contributions have been studied in a previous paper 
[2]. Therefore in the present paper we shall give our main attention to the interference contribution. 
Substituting (|2l]) into (|2^) and using ( p~5| ) we obtain for the single scattering contribution 

= TTT^Z I drdr'B(r r')A (r)A*(r') e "*W-*") 



167r 2 i? 2 e . 

[SgiSsj + n^jnl - 8 zi n j n z - 8-, j n i n 7 \ (22) 
Solving fl|), one obtains for the background potential 

M$ = -—%^- (23) 

Substituting (^) into (|2^) after some straightforward calculations we find for the single scattering 
contribution to the radiation intensity / = ^-Iu(R) [2], 

e 2 L z B(\k Q - k cos 9\) sin 2 9 u 2 
1 W = ^77— 7^— „ nl2 7X3 (2 4 ) 



2c [ 7 -2 + s in 2 6»fc 2 /fco] 2 ^ 
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— 1 /2 

where 7 = (1 — ev 2 /c 2 ) is the Lorentz factor of the relativistic particle in the medium, fco = 7, 
n z = cos 6 and L z is the system size in the z direction. As seen from ( p4| ) at ak <C 1 (B ~ const) , the 
forward and backward intensities are equal. When afc> 1, for relativistic particles &o — > k, 7 3> 1 the 
forward intensity (0 rs 0) is significantly larger than the backward intensity (0 ~ ir) because of the factor 
B. The diffusion contribution can be represented in the form [2] 

fc 2 



'^ B "> KiTT-'/?'^- 



dfdf'B(r - r')A (r)Ao(r') / drxdr^dr^dr^ 



e- ikft ^-^P(f 1 ,f 2 ,f 3 ,u)G(f f 3, r)G*{r', r 4 ) 
[fe^ij + riiUjul - 8 zi n j n z - S zj nin z ] 

where P is the diffusion propagator 



(25) 



n 



P(fi,f 2 ,r 3 ,r 4 ) = ^ 



(26) 



'"2 



Here the solid line denotes the averaged Green's function and the dotted one denotes the correlation 
function of random field. Respectively, the interference contribution has the form 

I%{R) = 16 ^2 R 2 £ J drdr'B(r-r')A (r)A*(r') J dndf 2 dr 3 dr A 
e-* kR ^- r ^P c (r u r 2 ,r 3 , r 4 )G(r 3 , f)G* (r ', r 4 ) 

[SziSzj + niUj-nl - 8 zi n j n z - 5 zj nin z ] (27) 



where 



ri T3 



(28) 



r 2 r A 



As follows from (jS^), due to time-reversal invariance (see, for example [7]) P is related to the diffusion 
propagator in the following manner 

P C (ri, r 2 , ri, ri) = P(n,rl, r^r^) (29) 

Thus the interference contribution as well as the diffusion one is determined by the diffusion propagator 
P. We now turn to finding it. Our consideration is similar to the three dimensional case [3]. It follows 
from (^6|) that P can be represented in the form 
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P(fi, r 2 , r 3 , r 4 ) = - r 2 )B(f 3 - r 4 )P(R', r*i - r 2 , r 3 - r 4 ) 



(30) 



where i?' = ^(js + — f\ — T2) and P satisfies the equation 



dp 
(2^ 



dq 



■f(q,K)B(p-q) 



P(K,p,q') = f(q',K) 



where 



M K) = G(q + K/2)G*(q- K/2) 



As usual, one needs to know P at K — > 0. In this limit the diffusion propagator has the form [3] 

p(k - o,p, = ^y^) 



where 



ImT,(q) 



{qK) 2 ImG(q) dq 



Im 2 S(g) (2tt) 3 

Substituting (|l^) and (14) into ( |34|) and calculating the integral in the limit fca > 1, we obtain 

20tt 1 



A{K) 



k 3K^l 2 +K 2 n P 



(31) 



(32) 



(33) 



(34) 



(35) 



In our previous paper [2] we investigated the special case K p = 0, which was sufficient for studying the 
diffusion contribution. Note that although we considered the case b» 1 all results are qualitatively 
correct also in the in the general case. When we know the form of diffusion propagator we can investigate 
the investigate the diffusion and interference contributions. First we consider the diffusion contribution. 
Substituting @ into (§|) and using @, (§|) and @, we obtain [2] 



I u (n) 



— (1 - n 2 z )A(K)k 2 lmT,(kn)L z x 



dq p 



B(\k + Jk 2 -q 2 \)+B[\k -Jk 2 -q 2 



(2tt)2 (qj + k 2 - k 2 f 



B{0) + B(2Jk 



(36) 



As seen from (|36|) , the main contribution to the integral over q p for relativistic particles 7 > 1, fco —> k 
give the values q p ~ 0. Accounting for this circumstance and using (|3^) we have from ( |36"| ) 

5 e 2 7 2 /' I 



6 ec 



3 • 2 a 

sin 6/ 



cost 



(37) 



When obtaining (|37]) we substitute as usually l/K 2 at if — » by L 2 (there with we assume that 
L z < L x , L y ). A more precise approach would be to solve the appropriate Schwarzschild-Milne equation 
for the present problem. This could bring overall prefactors of order unity [8,4]. For our present purpose 
we shall not be interested in these effects. It follows from ( |37| ) and (p4|), that I D /I° ~ (L z /l) 2 1. 
Thus when k\ cos9\l 3> 1 and I <C £ z the diffusion contribution is superior to the single scattering one. 
As one could expect, the forward and backward intensities in the diffusion contribution are equal to each 
other. Finally, we note the strong dependence of the diffusion contribution on the particle energy. 
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V. INTERFERENCE CONTRIBUTION 



Using ( J29| ) and Q30|) and changing the variables of integration by formulae 

xi = rl - fl, x 2 =fa — f^R' = -(r~5 + f 2 - f[ - fl),fl = fl (38) 

we find from (]2^ ) 

I C (H) = ^ ~ Uz } Ck [ dfdf"B(r-r')Ao(f)A*{f') [ dx^dR'dn 

32lT £ J J 

^ A -^^P{R',S 1 ,x 2 )B{S 1 )B{S 2 )G{R' + Tl±El+fl - f)G*(f' - r 4 ) (39) 

In the Fourier representation one finds from ( p9| ) 

T C(=i\ {l-n 2 z )ck 2 [ dqdqidq 2 dK -. 2 

7 (n) = 32.% J (27r )i2 g(g)|^(-*n-g-g)| 

B(gi)B(<f 2 )P(i?, fc, n + y - fcn + |- - 9 2)|G(fcn + /?)| 2 (40) 



Substituting (13), (|23|), ( |33| ) and (|3J) into (40), taking into account that the main contribution in the 
integral over K give the values K — » and sequentially integrating (^) using the Ward identity, we find 
from (Eol) 



jc _ 107re 2 L 2 |n z |(l - n z ) 2 B{\kcos8 + fc |) f dJT 1 (4 



£Cl J ^ (3K 2 Z + K 2 ) [(K p + kn p ) 2 + fc V 2 

Note that at ak > 1, B(2k)/B(0) ~ l/fc 2 a 2 < 1. Therefore, as follows from 0) in the Cherenkov limit 
the backward intensity (8 w 7r) is significantly larger than the forward intensity (8 « 0). This is the 
main characteristic feature of the interference contribution. It is analogous to the light back- scattering 
peak which occurs in propagation of light in the randomly inhomogeneous media [6] . Considering angles 
sin 8 ^> A/2nL p and calculating the integral over K in ( [fl] ) we have 

c SV^arctanv^ e 2 L z B(\kcos9 + k \)sin 2 9\ cos8\ 

' \") = o (42) 

^e c I 2 [k 2 sin 2 9 + k 2 j- 2 ] 

When obtaining (^2|) we cut the integral over K on the upper limit at 1/1. Note that the ratio A/2irL p , in 
the optical region, is of order ~ 10~ 4 — 10~ 5 . Therefore one can believe that the condition sin# ^> A/2irL p 
is always satisfied in the optical region. Comparing ( f42| ) with the single scattering contribution (]24|) we 
see that I c /I ~ l/(fc^) 2 <C 1. Although the interference contribution is small however it has quite 
different angular dependence. Accounting for the form of correlation function B (|l6|), it follows from 
( p!2| ) that the maximum of the radiation intensity for relativistic particles 7 > 1, k — > k, lies in the 
region of angles 9 ~ tt — 7 _1 . These angles are very close to the backward direction. On the other hand, 
the maximum of the single scattering contribution lies in the strongly forward range of angles 9 ~ 7 _1 . 
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VI. SUMMARY 



We have considered the influence of interference effects on the radiation of a charged particle passing 
through a stack of randomly spaced plates. It appears that interference contribution to the radiation 
intensity has a peak in the backward to particle motion direction. Though its value is small compared to 
the single scattering and diffusion contributions, it can be investigated experimentally. This is possible 
due its specific angular dependence. 
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